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We show that the Big Bang Observer (BBO), a proposed space-based gravitational-wave (GW) detector, 
would provide ultra-precise measurements of cosmological parameters. By detecting ~ 3 x 10"" compact-star 
binaries, and utilizing them as standard sirens, BBO would determine the Hubble constant to ~ 0.1%, and 
the dark energy parameters wo and Wa to ~ 0.01 and ~ 0.1, respectively. BBO's dark-energy figure-of-merit 
would be approximately an order of magnitude better than all other proposed, dedicated dark energy missions. 

To date, BBO has been designed with the primary goal of searching for gravitational waves from inflation, 
down to the level Qaw ~ 10~^^; this requirement determines BBO's frequency band (deci-Hz) and its sensi- 
tivity requirement (strain measured to ~ 10"^''). To observe an inflationary GW background, BBO would first 
have to detect and subtract out ~ 3 x IC merging compact-star binaries, out to a redshift z ~ 5. It is precisely 
this carefully measured "foreground" which would enable high-precision cosmology. BBO would determine the 
luminosity distance to each binary to ^ percent accuracy. In addition, BBO's angular resolution would be suf- 
ficient to uniquely identify the host galaxy for the majority of binaries; a coordinated optical/infrared observing 
campaign could obtain the redshifts. Combining the GW-derived distances and the electromagnetically-derived 
redshifts for such a large sample of objects, out to such high redshift, naturally leads to extraordinarily tight 
constraints on cosmological parameters. We emphasize that such "standard siren" measurements of cosmology 
avoid many of the systematic errors associated with other techniques: GWs offer a physics-based, absolute 
measurement of distance. In addition, we show that BBO would also serve as an exceptionally powerful grav- 
itational lensing mission, and we briefly discuss other astronomical uses of BBO, including providing an early 
warning system for all short/hard gamma-ray bursts. 

PACS numbers: 04.30-w,04.80.Nn,95.36.+x,95.35.+d,98.62.Sb,98.80.Es 



I. INTRODUCTION 

Improving our understanding of the dark energy respon- 
sible for the observed accelerating expansion of the Uni- 
verse is one of the foremost challenges in physics. Our cur- 
rent theoretical models are exceedingly inadequate. In ad- 
dition, given its extremely low energy density, there are no 
plausible scenarios for the direct detection of dark energy. 
Thus progress on this critical issue must be made through 
indirect observations. In this paper we show that the Big 
Bang Observer (BBO), a proposed space-based gravitational- 
wave (GW) mission designed primarily to search for inflation- 
generated stochastic GWs in the band 0.03Hz-3Hz [1], 
would also be an ultra-precise cosmology mission, measuring 
the Hubble constant Hq and the dark energy parameters wq 
and Wa far more accurately than other proposed dark energy 
missions. 

BBO has been proposed as a follow-on mission to the 
Laser Interferometer Space Antenna (LISA), which will be 
most sensitive to GWs in the band ^ 10^'*-^-10~^-^Hz. In 
the LISA band, any primordial GWs from standard infla- 
tion are likely to be buried under the GW foreground from 
the short-period white dwarf-white dwarf (WD-WD) binaries 
in the universe. The WD-WD foreground is practically ab- 
sent in the BBO band, since the WD-WD contribution falls 
very rapidly for / > 0.01 Hz, and disappears entirely for 
/ > 0.25 Hz, where the most massive (and hence smallest) 
WDs merge lH] El O. Instead, the dominant astrophysical 
foreground in the BBO band is mergers of compact bina- 
ries composed of neutron stars (NSs) or black holes (BHs); 



i.e., NS-NS, NS-BH, and BH-BH binaries. BBO would be 
sufficiently sensitive to individually detect and subtract out 
essentially every merging compact binaries out to high red- 
shift, thereby uncovering any primordial GW background in 
its band that has energy density flaw ^ 10"^'' HHHIll. 

As far as the search for primordial GWs is concerned, merg- 
ing compact binaries represent a foreground that must first 
be removed. In this paper we show that this foreground is a 
cosmological gold mine, allowing astronomers to measure the 
expansion history of the Universe (out to at least z ~ 5) far 
more accurately than the most ambitious currently-proposed 
cosmology missions. The basic argument is this: GW de- 
tectors in general, and LISA and BBO in particular, will 
provide high-accuracy measurements of the luminosity dis- 
tances, D]^, to detected merging binaries. However, the GW 
signals contain (almost) no information regarding the source 
redshifts. Thus the situation in GW astronomy is the exact 
reverse of optical/electromagnetic astronomy: accurate dis- 
tances will be relatively easy to come by, while determining 
redshifts will be much more challenging. For both the ground- 
based LIGO/VIRGO network and the space-based LISA, the 
angular resolution for detected binaries will typically be sev- 
eral degrees Id, J], so the error box on the sky would contain 
^ ^q5-6 possible host galaxies per event. Therefore plans for 
doing cosmology with GW observations have usually hinged 
on finding some electromagnetic outburst associated with the 
GW events, in order to identify the host galaxy and obtain 
its redshift [8l (though see Finn and Chernoff f9l, MacLeod 
and Hogan |10|, and Seto et al. fTTl for suggestions on how 
to evade this requirement). GW-detected merging binaries for 
which one can also determine a redshift have been dubbed 



2 



"gold-plated" binaries; Holz and Hughes have shown that, 
with LISA, even a handful of such "gold-plated" detections 
of massive black hole binaries could make significant contri- 
butions to cosmology |12|. 

BBO should detect ~ 10^ merging NS-NS binaries per 
year ||4l, out to z '--^ 5. BBO's angular resolution for NS-NS 
mergers will typically be a few arcsec: a small enough error 
box to uniquely identify the host galaxy in most cases. This is 
especially true since one can also use the source's luminosity 
distance, Dl (measured to several percent), and even a crude 
Dl-z relation, to rule out galaxies at the right position on the 
sky but with significantly different redshifts (e.g., differing by 
> 10%). While a unique host galaxy may not be identifiable 
in very dense galaxy clusters, in such cases one can substitute 
the average redshift of the cluster; this will lead to a typical 
redshift error of order Az ~ 2 x 10^'^, which is negligible for 
our purposes. Thus — in stark contrast with LISA and ground- 
based GW detectors — BBO will detect ^ 10^ "gold-plated" 
binaries per year! BBO will measure the luminosity distance 
to each NS-NS binary with a relative accuracy of several per- 
cent. For example, for a NS-NS binary at z = 1.5, the median 
distance error due to detector noise will be ~ 2%, while the 
distance error due to weak lensing (WL) will be ~ 7%. For 
a large sample of sources, both of these errors "average out", 
and the extent to which they do not average out can be read- 
ily modeled and accounted for By contrast, the ambitious, 
space-based dark-energy mission SNAP would be expected to 
observe roughly 2,000 SNe over the lifetime of the mission, 
out to a maximum redshift of z ~ 1.7 ifTSI . BBO provides 
an overwhelmingly larger and deeper data set, with each in- 
dividual distance measured significantly more precisely. It is 
thus to be expected that BBO will vastly outperform other pro- 
posed dark-energy missions; in fact, BBO might well provide 
better constraints than all other proposed dark energy mis- 
sions combined. As we show below, with the inclusion of 
an (expected) Planck prior, BBO data should measure Hq to 
- 0.1%, Wo to 0.01 and Wa to 0.1. 

Perhaps even more importantly, we argue that the system- 
atic errors associated with GW detections are generally much 
smaller, and much easier to characterize, than with any other 
proposed methods (e.g., weak lensing, supemovae, or baryon 
acoustic oscillations). Type la supernovae (SNe) are arguably 
the state of the art in cosmological distance measurement. 
The intrinsic luminosity of a type la SN can be empirically 
calibrated to roughly 10%. The physics underlying this cali- 
bration is only poorly understood, and possible evolutionary 
systematic effects are a grave cause for concern (see, e.g., 
lfT4l [TSl [T6ll ). By contrast, the compact binaries detectable 
by BBO are exceptionally simple sources. At the orbital sep- 
arations at which BBO observes them, the compact objects 
can be treated as point masses (with spin), whose dynamics 
are very accurately described by the post-Newtonian approx- 
imation. Systematic distance errors arising from the detector 
itself will also be negligible, since BBO (like LISA) is fun- 
damentally self-calibrating 1 17|. The optical scheme is some- 
what complicated in practice, but in essence one is measuring 
the spacecraft arm lengths (more precisely, small time-varying 
changes in the differences between spacecraft arm lengths) in 



units of the laser wavelength. For LISA, the laser wavelength 
(or equivalently its frequency) will be known to an accuracy 
of ^ 10^^-10^^; this small uncertainty is expected to dom- 
inate LISA'S calibration error, which is therefore also at the 
^ 10^^-10^^ level. (For the same reason, LISA Pathfinder's 
calibration is expected to be accurate to better than 10^^ IIS].) 
One could reduce this small inaccuracy by, say, stabilizing the 
laser to a transition line in an iodine gas cell [ 19 1, if there were 
sufficient motivation. However there is one caveat concerning 
BBO's calibration accuracy, which we discuss in Sec. |IIIC| 

In addition to determining the equation of state of dark en- 
ergy, BBO should also be an excellent probe of the growth of 
structure. BBO's measurements of gravitational lensing will 
be as sensitive as dedicated weak lensing missions. The ba- 
sic idea is that once a Dj^-z relation has been derived from 
the entire NS-NS sample, the dispersion of the ^ 3 x 10^ 
NS-NS data points about this curve is dominated by magni- 
fication by gravitational lensing. That is, one has obtained 
~ 3 X 10^ independent measurements of the lensing mag- 
nification along different lines-of-sight, one for each binary. 
We show below that the typical SNR for each measurement is 

3.4, and theref ore the S NR for the whole NS-NS popula- 
tion is - 3.4 X V3 X 105 - 2 X 10^. Although the NS-NS 
merger rate most likely exceeds the BH-BH rate, the lensing 
SNR from the BH-BH population may exceed that of NS-NS 
mergers. The BH-BH merger rate is poorly constrained by 
observation; a reasonable estimate from population synthe- 
sis models is that the BH-BH rate is a factor ^ 20 smaller 
than the NS-NS rate |20|. However, the typical SNR for each 
BH-BH merger is larger by a factor ~ 5.3 (for our fiducial 
NS and BH masses). We thus estimate a lensing SNR for the 
whole BH-BH population of - 2.2 x 10^, and a total SNR for 
both populations combined of 3 x 10'^. This total lensing 
SNR could be significantly higher, if the BH-BH merger rate 
is near the high end of the estimated range. (Of course, ad- 
vanced ground-based GW detectors will measure these rates, 
for z < 0.4, many years before BBO flies.) In Sec. Ill B we in- 



vestigate BBO's sensitivity as a gravitational lensing mission, 
and compare it with other lensing missions. 

While obtaining redshifts for 3 x 10^ host galaxies would 
be a highly ambitious goal at present, we are optimistic that it 
will be far less daunting by the time BBO flies. By then LSST 
may already have determined photometric redshifts (accurate 
to ^ 2-3%) for a large fraction of the host galaxies in ^ 1/3 
of the sky. In addition, there are many proposed wide-field 
spectroscopic surveys; for example, BigBOSS would mea- 
sure ^ 5 million spectroscopic redshifts/yr (4,000 at a time) 
for galaxies in the range 0.2 < z < 3.5, over an area of 
14,000deg2 [21 1. The success of BBO as a cosmological 
probe is dependent upon the determination of host redshifts. 
It will thus be critical to secure the necessary optical resources 
and develop an efficient strategy for obtaining redshifts for a 
large sample of binary host galaxies. 

After this work was mostly completed we found some brief 
remarks in the literature that partly anticipate our results. 
Crowder & Cornish |22 | point out (in one sentence) that BBO 
should be able to localize the host galaxy for most observed 
compact-binary mergers, but the profound implications of this 



fact for physical cosmology are not developed. One slide in 
a workshop talk on Decigo by Seto [23] lists in brief bullets 
that Decigo measurements of NS-NS mergers could be used 
to probe dark energy though the D^-z relation, and that some 
of these detections will likely be GW-counterparts to gamma- 
ray bursts, allowing one to get a redshift. But these slides give 
no estimates of the resulting cosmological constraints, nor 
any sense of the potentially revolutionary implications of this 
observation (perhaps because the current version of Decigo 
would likely have relatively poor calibration accuracy, and so 
the cosmological measurements would suffer from large sys- 
tematics; see Sec. IIICi. Seto et al. |11| also proposed us- 



ing deci-Hz GW detections of inspkalling NS-NS binaries 
to measure the universe's expansion, using a fundamentally 
different approach: observing the small fractional change in 
source redshift over the course of the inspiral (detected as a 
small time-varying phase shift). Because this phase shift is 
such a small effect, the cosmological constraints that these 
authors estimated as obtainable are orders of magnitude less 
sensitive than the constraints we find (for comparable mission 
sensitivities). 




FIG. 1: Big-Bang Observer (BBO) consists of four LISA-like trian- 
gular constellations orbiting the Sun at 1 AU. The GW background 
is measured by cross-correlating the outputs of the two overlapping 
constellations, while time-of-flight across the Solar System gives 
BBO its angular resolution. A schematic of Decigo (a Japanese 
proposal similar to BBO) would be almost identical, except that the 
constellations are 50 times smaller than BBO's, and their arms form 
Fabry-Perot cavities. 



II. OVERVIEW OF BBO, DECIGO, AND THE NS-NS 
MERGER FOREGROUND 

A. BBO 



In this paper we focus on BBO, although it is to be noted 
that the Japanese GW community is proposing a similar mis- 
sion called Decigo. Our basic conclusions apply to any deci- 
Hz GW mission of roughly BBO-level sensitivity, including a 
Decigo-like mission, so long as the design ensures excellent 
calibration accuracy. We discuss this further in Sec. |II A| 



The rest of this paper is organized as follows. In Sec. II 
we give a brief overview of the BBO mission, its design 
sensitivity, the foreground produced by NS-NS binaries, and 
how accurately BBO can measure the distance and sky lo- 
cation of inspiralling compact binaries. In Sec. Ill we de- 
rive the cosmological measurement accuracy obtainable by 
BBO from 3 x 10^ NS-NS binaries. We also investi- 
gate BBO's performance as a weak lensing mission, and dis- 
cuss some caveats accompanying our conclusions. In Section 
IV we briefly discuss several other astrophysical uses of the 
BBO compact-binary data: as an early warning system for 
all short/hard gamma-ray bursts, in searches for the earliest 
intermediate-mass black hole mergers, and in studies of hun- 
dreds of strongly lensed GW sources. Our conclusions and 
plans for future work are outlined in Sec. V. 



We use units in which G = c = 1; everything can be mea- 
sured in the fundamental unit of seconds. However, for the 
sake of familiarity, we also sometimes express quantities in 
terms of yr, Mpc, or AIq, which are related to our fundamen- 
tal unit by 1 yr = 3.1556 x lO'^ s, 1 Mpc = 1.029 x 10" s, 
and IMq — 4.926 x 10^® s. For concreteness, in our simu- 
lations we adopt the following fiducial values for the cosmo- 
logical parameters: = 70kms^^ Mpc^^, = 0.3, and 
= 0.7. 



BBO will be most sensitive in the band ^ 0.03-3 Hz, which 
is dictated by BBO's main design goal: to detect primordial 
GWs generated by inflation. BBO has been proposed as a 
follow-on mission to the Laser Interferometer Space Antenna, 
which will operate at - lO^'^-^-lO^^-^Hz. In the BBO band, 
the dominant astrophysical foreground will be mergers of NS- 
NS, NS-BH, and BH-BH binaries. The BBO mission is de- 
signed so that essentially each and every merging compact bi- 
naries in the observable universe (i.e., on BBO's past light 
cone) can be detected and subtracted out. Cutler & Harms 0| 
have presented a detection/subtraction algorithm and showed 
by an analytical calculation that it should work; Harms et 
al. |5| implemented the algorithm and demonstrated that it 
worked well on simulated BBO data, albeit thek demonstra- 
tion was on a much reduced data set because of computational 
limitations. A very different algorithm for detecting and sub- 
tracting the binary signals from BBO data is also being devel- 
oped by N. Kanda and collaborators (unpublished). 

The current BBO design calls for four constellations of 
three satellites each, all following heliocentric orbits at a dis- 
tance of 1 AU from the Sun (see Fig. [T]). Each 3-satellite 
constellation can be thought of as a highly sensitive "mini- 
LISA", since the armlengths for each constellation are two 
orders of magnitude smaller than LISA'S, and BBO's sensi- 
tivity band is correspondingly two orders of magnitude higher 
in frequency. Two of the constellations overlap to form a "star 
of David"; the other two are ahead and behind by 27r/3 ra- 
dians, respectively. Briefly, the idea behind this orbital ge- 
ometry is that the energy density of the primordial GW back- 
ground, riGw(/), will be measured by cross-correlating the 
outputs of the two overlapping constellations in the star of 
David (much as LIGO attempts to measure ilcwif) by cross- 
coiTelating the outputs of the Livingston and Hanford inter- 
ferometers |25|). The other two constellations give BBO its 
angular resolution, which is useful for characterizing and re- 
moving the merging compact binary foreground. The source's 
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angular position on the sky is mostly determined by triangu- 
lation, using the differences in arrival times of the GWs at the 
different constellations. 

As explained in Cutler & Harms |4| , for compact-binary 
mergers the science output of each mini-LISA is in practice 
equivalent to the output of two synthetic Michelson detectors, 
represented by the time-delay interferometry (TDI) variables 
X and {Y — Z)/\/3. We can therefore regard BBO, which 
is made up of 4 mini-LISAs, as formally equivalent to 8 syn- 
thetic Michelson interferometers. To construct the instrumen- 
tal noise curve, 5h(/), of each of these synthetic Michelsons, 
we use Larson's on-line "Sensitivity curve generator" |26| and 
plug in BBO's instrumental parameters, which are taken from 
the BBO Concept Study 1 1 1 and also listed in Table [l] These 
parameters will be subject to change as the mission evolves, 
but for now they provide a convenient baseline. The BBO 
Concept Study lUl also lists parameters for less and more am- 
bitious versions of the BBO mission, referred to as "BBO-lite" 
and "BBO-grand", respectively, but in this paper we restrict 
attention to the intermediate version, or "standard BBO". In 
using the on-line generator, we have specified that the high- 
frequency part of S'h is 4 times larger than the contribution 
from photon shot noise alone; this factor 4 accounts for high- 
frequency noise components other than shot noise, such as 
beam pointing jitter and stray light effects. This is the same 
choice made in Fig. 1 of the BBO proposal 1 1 1, and is consis- 
tent with the standard assumptions made in drawing the LISA 
noise curve. This BBO instrumental noise curve is shown in 
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FIG. 2: Amplitude of BBO's instrumental noise, \/ fSh{f ), com- 
pared to the amplitude of the (pre-subtraction) NS binary foreground 
(plotted for ho — 10^^ Mpc^^yr^^) and the sought-for cosmic GW 
background (plotted for Qcwif) ~ 10^^''). To reveal this cosmic 
background, the NS foreground must be subtracted off, with frac- 
tional residual of < 10"^'^. 



suing a BBO-style mission, but with a somewhat less ambi- 
tious sensitivity goal, might be advisable from a cost/benefit 
standpoint; Decigo may have a comparable cosmological 
reach to BBO, if designed to ensure excellent calibration ac- 
curacy. In follow-up work, we plan to investigate how the 
science payoff from a deci-Hz GW mission varies with its sen- 
sitivity. 



BBO is seen as a follow-on mission to LISA in the U.S. 
and Europe, but in the Japanese GW community there is a 
strong push to launch a deci-Hz GW mission first. The current 
plan is for Decigo to be a factor ~ 2-3 less sensitive than 
BBO, but for it to launch earlier, in ~ 2024. A small Decigo 
precursor mission, Decigo Pathfinder, is among the final two 
missions competing for a launch slot in ^ 2012; JAXA (the 
Japanese space agency) is scheduled to decide in mid-2009 
which mission gets this slot. 

Our research to date has concentrated on BBO, but it would 
be straightforward to generalize our BBO analysis to a mis- 
sion with Decigo-level sensitivity. Indeed, we expect that pur- 





Symbol 


Value 


Laser power 


P 


300 W 


Mirror diameter 


D 


3.5m 


Optical efficiency 


e 


0.3 


Arm length 


L 


5 ■ lO^m 


Wavelength of laser light 


X 


0.5 /im 


Acceleration noise 


\^ Sa,cc 


3- lQ-^^m/(s2VHz) 



B. The NS-NS merger rate over time 

BBO would be able to observe BH-BH and BH-NS merg- 
ers to significantly higher redshifts than NS-NS mergers, but 
the rates for BH-BH and BH-NS mergers are more uncertain 
(and probably a factor ^ 20 lower) than NS-NS mergers, so 
our discussion in this section will focus mostly on the NS-NS 
case. The extension of our work to NS-BH and BH-BH bina- 
ries is straightforward. 

We denote the NS-NS merger rate (per unit proper time, per 
unit co-moving volume) at redshift z by n{z). It is convenient 
to regard n(z) as the product of two factors: 



h{z) = fiQ ■ r(z) 



(1) 



where ho is the merger rate today and r{z) encapsulates the 
rate's time-evolution. For r{z), we adopt the following piece- 
wise linear fit to the rate evolution estimated in ETl : 



r{z) 



l + 2z 

1(5 -^) 




z < 1 

1 < z < 5 

z > 5 



(2) 



TABLE I: BBO parameters 



The current NS-NS merger rate, no, is also usefully regarded 
as the product of two factors: the current merger rate in the 
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Redshift z 

FIG. 3: The total number of NS-NS mergers closer than redshift z. 
The results are normalized to a 3-yr observation period and no = 

10"''Mpc"^yr-\ 



Milky Way, and a factor that extrapolates from the Milky Way 
rate to the average rate in the universe. The NS-NS merger 
rate in the Milky Way has been estimated by several authors; 
it is still highly uncertain, but most estimates are in the range 
10-M0"'*yr-i EHllISl'SOl. To extrapolate to the rest of the 
universe, Kalogera et al. |29 1 estimate that one should multi- 
ply the Milky Way rate by 1.1-1.6 x 10"^ • /170 Mpc"^. This 
factor is obtained by extrapolating from the B-band luminos- 
ity density of the universe, and it is only a little larger than the 
extrapolation factor derived by Phinney in OTI . For the esti- 
mates in this paper we adopt a rough geometric mean of these 
rates: no = 10~^Mpc~^yr~^. 

How many NS-NS merger events, AiVm, enter the BBO 
band during some observation time, Atq? Integrating the con- 
tributions from all redshifts, the rate is given by ||4]: 

AiV„, = 3.0 . 10^ ( ^\ ( , 3 ) . (3) 

Based on our fiducial merger rate Fig.jsjplots the num- 
ber of observable mergers during three years of observation 
that occur closer than (any given) redshift z. We see that 
roughly 15% are at z < 1, the median redshift is z ^ 1.6, 
and roughly two-thirds are between z — \ and 3. This distri- 
bution is well suited for probing the evolution of dark energy, 
since it fully samples the evolution history during the transi- 
tion from the matter dominated to the dark energy dominated 
era |32|. 

The time required for a NS-NS inspiral signal to sweep 
through the BBO band will typically be comparable to BBO's 
lifetime. More specifically, the time remaining until merger, 
from the moment the GW frequency sweeps through a given 
frequency, /, is (to lowest post-Newtonian order) 

\ 1.22M0 j viHzy ^ ' 

where M. = /.t'^/^Af^/^ is the so-called "chirp mass" of the bi- 
nary (with M the binary's total mass and /i its reduced mass). 



For two 1.4MoNSs,/« 0.205 Hz, 0.136 Hz, and 0.112 Hz 
at one year, three years, and five years before merger, respec- 
tively. 

C. BBO's accuracy of parameter estimation 

Although it is straightforward to calculate BBO's angu- 
lar resolution and distance measurement accuracy using the 
Fisher matrix approximation, for pedagogical reasons we be- 
gin with some simple, back-of-the-envelope estimates. If the 
source's luminosity distance, D^, were the only unknown pa- 
rameter, it could be measured to a relative accuracy of 1/SNR, 
where SNR is the (amplitude) signal-to-noise ratio of the de- 
tection. For a merging NS-NS binary at z = 1, the median 
SNR will be roughly 180, suggesting I^Dl/Dl ~ 0.6%. 
This is a lower limit on the eiTor, since correlations with 
other unknown parameters increase the uncertainty. Experi- 
ence from similar problems in GW data analysis suggests an 
increase by a factor of ^ 2, or ^ 1% distance uncertainty at 
z=l. 

The source's angular position on the sky will be inferred 
mostly by triangulation, based on the GW time of flight be- 
tween the different mini-LISA constellations. This idealiza- 
tion suggests the estimate (l/SNR)(l/27r500/o) 
1 arcscc, where /o ~ 0.3 Hz is BBO's most sensitive fre- 
quency and we have used 1 AU « 500 s in geometric units. 
This estimate neglects both the information about the source's 
sky location that is encoded in the waveform by the time- 
varying antenna patterns of the constellations, as well as the 
increased uncertainty that results from correlations with the 
other unknown parameters. 

These back-of-the-envelope values turn out to be reason- 
ably good estimates of BBO's distance and angular position 
accuracies. We now turn to a more careful. Fisher matrix cal- 
culation of BBO's parameter estimation accuracy. The ap- 
plication of Fisher matrix techniques to LISA measurements 
of inspiralling binaries was demonstrated in detail in Cut- 
ler f7|. The generalization to BBO, which is just four mini- 
LISAs with higher sensitivity, is straightforward; our exposi- 
tion will therefore be biief, and we refer the reader to |7 1 for 
more detail. As mentioned previously, we can regard BBO's 
output as formally equivalent to that of 8 independent, syn- 
thetic Michelson interferometers; we represent BBO's output 
as Sa{t), for a = 1 • • • 8. We use s to abstractly represent 
these 8 time-series. For simplicity we assume that the de- 
tector noise is stationary, Gaussian, and the same for all four 
mini-LISA's (in practice, we expect that the noise levels will 
be somewhat different and slowly time-varying, in a manner 
that we can fit for). Under these assumptions, we obtain the 
following natural inner product on the vector space of signals. 
Given two signals g and k we define (g | k) by 

where ga{f) and ka{f) are the Fourier transforms of ga{t) 
and ka{t), respectively. We follow the usual convention of 
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taking Sh{f) to be the single-sided, sky-averaged noise spec- 
trum for each synthetic Michelson interferometer. The fac- 
tor 3/20 in Eq. (jsj) is the product of a factor 1/5 due to the 
sky-average convention and a factor 3/4 = sin^(7r/3) arising 
from the 7r/3 angle between the arms in each constellation; 
see Sec. V. A of Barack & Cutler ll33l for a fuller explanation. 
In this notation, the rms SNR for any waveform h is 

SNR[h] = (h|h)i/2. (6) 

For a given incident gravitational wave, different realizations 
of the noise will give rise to somewhat different best-fit pa- 
rameters. However, for large SNR, the best-fit parameters will 
have a Gaussian distribution centered on the correct values. 
Specifically, let be the "true" values of the physical pa- 
rameters, and let + AA^ be the best fit parameters in the 
presence of some realization of the noise. Then for large SNR, 
the parameter-estimation errors AA*^ have a nearly Gaussian 
probability distribution whose covariance matrix is given by 
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AAA'AA" = (r-i)'"" ( 1 + ©(SNR) 



(7) 



FIG. 4: BBO's angular resolution as a function of redshift, z. The 
three curves show BBO's median la angular resolution for three 
fiducial types of merging compact binaries: BH-BH, BH-NS, and 
NS-NS. 



where the overline " — 
where F^i^ is the Fisher matrix, defined by 



means "expectation value", and j^e waveform phase is given by 



r = 
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Cutler and Harms 0] have shown that the effects of orbital 
eccentricity on the NS-NS GW signal will typically amount 
to less than one radian of phase over the entire ^10^ radians 
of observed inspiral, while the precession of the orbital plane 
due to the Lense-Thirring effect will typically be < 10^"^ ra- 
dians in the BBO band. We can therefore model the binaries 
as quasi-circular, and neglect spin-precession effects (but we 
do include spin-orbit effects on the waveform phase). Our sig- 
nal waveform, ha{t), thus depends on 10 physical parameters 
parameters describing the binary: Mi, M2, (3, Og, 4>s, 4'l, 
(f)c, tc, and D^. Here Mi and M2 are the "redshifted mass" 
(Mi = (1 + z)Mi), /3 is a spin-orbit coupling parameter de- 
fined in Cutler & Flanagan |6|, (6*5, (pg) give the direction 
to the source, {9l, describe the orientation of the orbital 
plane, Dl is the luminosity distance to the source, tc is the 
time of merger, and 0c is a constant of integration in the evo- 
lution of the binary's orbital phase. Our signal model is given 
by: 

hM) ^ ^^rV6A4i)e<*(/)-<W--= W) (/ > 0) 

(9) 
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where the PN expansion parameter a;(/) is defined by 
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x{!) = ( ^M(l + z)! ' 



-5/3 



,(12) 
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Equations for the modulation factor Aa{t{f)) and "polariza- 
tion phase" (faitif)) (both arising from the rotation of each 
mini-LISA constellation as it orbits the Sun) as well as for 
the "Doppler phase" <I>^ (due to time delay between the pas- 
sage of a particular wavefront over the Solar System barycen- 
ter and its passage over each constellation) are given explic- 
itly in Sees. Ill and V of Cutler |7|. The exact expressions 
depend on the two angles {(3q, ao) that describe each constel- 
lation's position around the Sun, and the orientation of each 
detector-triangle within its plane at some fiducial time, to. For 
definiteness, for the four mini-LISA's we choose: (Pq, a^) = 
(0, 0), (0, tt), (27r/3, 27r/3), and (47r/3, 47r/3). 

The uncertainty in the source's angular position, AUg (in 
solid angle), is given by | 



where M 
and 



Ml + M2, /i = M1M2/M, M 



Ails = 27r 



^{A^xsriA<Psr~{AfisAq^sr 



(14) 
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The relation between time and frequency, t = t{f), is given 
through 0([u/c] 3) by |6| 
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The 2tt factor on the right-hand side (RHS) of Eq. is con- 
ventional; with this definition, the probability that the source 
lies outside an (appropriately shaped) error ellipse enclosing 
solid angle An is g-'^^/'^^s. That is, as defined above, A^s 
is very good approximation to the size of the Icr error ellipse. 

BBO's median SNR angular resolution and distance accu- 
racy (both Ict) for NS and BH mergers and a range of z are 
shown in Figs.|4]and |5] respectively. These figures were pro- 
duced as follows. Each NS was taken to have mass 1.4Af0, 
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and each BH to have mass IOMq. For BBO's (sky-averaged) 
noise spectral density, Sh{f), we adopted the fitting function 



Shif) = e.lSxlO-^V"-* + 1.95x10--*^ + 1.2xlO-■*V^ 

(15) 

where / is in units of Hz. For each z, we chose 250 ran- 
dom angle sets (6s, (j>s, ^l, 0l), and computed the SNR, 
the Fisher matrix, and its inverse. Since the Fisher ma- 
trices are nearly degenerate, we tested robustness by using 
both Matlab's standard matrix inversion function and Matlab's 
Cholesky-factorization inversion routine; these were found to 
give essentially identical results. While we argued above that 
the spin-spin coupling will have a negligible impact on NS- 
NS waveforms for typical cases, as a further test of robustness 
we added an additional spin-spin parameter (usually denoted 
"cr" in the literature), and recomputed parameter estimation 
accuracies. The results shown in Figs. |4] and |5] turn out to be 
essentially independent of the presence/absence of a spin-spin 
term in the waveform model. Note that the results in Figs. |4] 
and|5]are in reasonably good agreement with the z — 3 results 
in Table V of Cornish & Crowder |22 | (considering that we 
model the high frequency part of BBO's noise curve some- 
what differently). 

Let us suppose that BBO identifies a binary system some- 
where in the universe. We now determine the number of po- 
tential host galaxies for the binary to be found in the BBO 
error volume. We closely follow the approach of Holz & 
Hughes |[T2l . updating their value for the projected number 
density to the Hubble Ultra Deep Field number: dN/dil — 
1, 000 galaxies/arcmin^ ll34ll . Since BBO measures distances 
at the percent level, the depth of the BBO error box is domi- 
nated by the distance uncertainty due to gravitational lensing. 
Following Eqs. 6-8 and Fig. 8 of lfT2l . we calculate the total 
number of galaxies in the BBO error box, per arcmin^. We 
then multiply this by the size of the BBO error box (shown in 
Fig. |4]), to arrive at the total number of galaxies in the BBO 
error box, as a function of redshift. This result is shown in 
Fig. |6] The largest number of galaxies in a BBO error box 
is in the case of NS-NS binaries at z ~ 1.5, but even in 
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FIG. 6: Number of galaxies in the BBO error cube, as a function of 
redshift. Even in the worst case, there is less than one galaxy within 
la of a given binary on the sky, and therefore it should be possible 
to robustly identify the unique host galaxy. 



this case there is less than "half" of a galaxy present. Thus 
even at the "worst" redshift, the median occupation fraction 
is less than one — it will be possible to identify the unique 
host for the majority of BBO sources, and hence associate 
the appropriate redshift for the majority of distance measure- 
ments. This is in contrast to the case of LIGO or LISA, where 
the error boxes are large enough that associated electromag- 
netic activity (such as a gamma-ray burst, or activity asso- 
ciated with a supermassive binary black hole merger) is re- 
quired to uniquely identify the counterpart lfT2ll35l[36ll . Obvi- 
ating the need for an independent identification of the counter- 
part sharply increases the expected number of usable standard 
sirens, and hence significantly improves the accuracy of the 
cosmological measurements. Galaxy misidentifications will 
generally be seen as large outliers, and thus their influence 
can be mitigated by the use of robust statistics, such as the 
Hough Transform. 
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FIG. 5: BBO's median la distance accuracy as a function of redshift, 
z, for merging BH-BH, BH-NS, and NS-NS binaries. 



III. ULTRA-HIGH PRECISION COSMOLOGICAL 
PARAMETERS FROM BBO 

A. The Dl-z relation 

We begin by considering BBO's measurements of the lu- 
minosity distance-redshift relation (see Fig. |7]i. This relation 
is a direct measure of the evolution history of the Universe: 
redshift provides the size of the Universe at emission, and lu- 
minosity distance provides the time since emission. Thus a 
precise measurement of this relation is sensitive to dark en- 
ergy; indeed, it is this method that enabled the initial discovery 
of the accelerating expansion of the Universe now associated 
with dark energy. 

We consider a fiducial population of 2.5 x 10^ NS/NS bi- 
naries distributed according to Eq. |2] out to z = 3. We as- 
sume that the distance measurement errors due to detector 
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noise for each individual binary are those shown in Fig. |4] 
Because BBO does such an exquisite measurement of dis- 
tance, the errors on the true distance to a given binary will 
be dominated by the effects of gravitational lensing magnifi- 
cation ll37l l38l . We incorporate the lensing errors following 
the approach of ll39l . which is entirely appropriate given the 
very high-number statistics we are considering. For each in- 
dividual binary we take the dispersion in flux due to lensing 
to be given by crionsing = 0.088z (see Eq. 9 of ll39l ). We have 
explicitly checked that this approach is equivalent to draw- 
ing magnification values from the full, non-Gaussian lensing 
probability distribution functions derived in |38|. We assume 
that the sky localization is sufficient for the identification of 
a unique host galaxy (and hence redshift) for each binary (as 
in Fig. |4|i. The redshift determination will need to be done 
independently of BBO, in the electromagnetic band. While 
in practice there will be some host galaxy misidentifications, 
for simplicity in this study we assume that perfect redshifts 
have been obtained for all of our sources. (This simplification 
is partly based on our belief that a robust cosmological pa- 
rameter estimation method will substantially mitigate the ef- 
fects of a fractionally small set of misidentifications — enough 
so that in estimating BBO's performance, to a first approx- 
imation it is reasonable to neglect them.) We Monte Carlo 
generate populations of observed binaries, and then for each 
population we determine the best-fit cosmological parameters 
(varying the number of free parameters of interest). We re- 
peat this procedure for a large (> 10^) number of runs, and 
plot the resulting error contours. In what follows, the la con- 
tours contain 68.3% of the best-fit values, and the 2a contours 
contain 95.5% of the models. 

We follow the common convention of parameterizing the 
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FIG. 7: Distance versus redshift for a sample BBO binary popu- 
lation. Distance is shown as distance modulus, and includes both 
BBO errors and gravitational lensing. The red curve is the true lumi- 
nosity distance-redshift relation. Notice that lensing causes a small 
number of binaries to become tremendously magnified (to lower 
distance modulus), but there is a lower limit to the amount of de- 
magnification. 



dark-energy equation of state in the two-parameter form BOl 

2; 

w{z) = Wo+Wa-IJ-^^. (16) 

We fit each data set to five cosmological parameters: the Hub- 
ble constant Hq = /i x 100 km/s/Mpc, the dark-matter density 
r^m, the dark-energy density ilr^, and the dark-energy phe- 
nomenological parameters wq and Wa- As is standard in as- 
sessing the power of proposed cosmology missions, we in- 
clude a forecasted Planck CMB prior, which constrains the 
angular diameter distance at z = 1080 to 0.01%, and con- 
strains flrnh'^ to 1% Bni42l . 

Fig.[8]shows the resulting constraints on h and ri,„, assum- 
ing our fiducial population of binaries, and a 5-parameter fit to 
the data. We find that BBO will measure the Hubble constant 
to ^ 0.1%, even when marginalizing over two dark-energy 
parameters. For comparison, the Hubble Key Project (one of 
the major goals of the Hubble Space Telescope) arrived at a 
value of 72 ± 8 km/sec/Mpc (> 10% error) for Hq |43|. It 
is to be noted that, if we fit the data to a ACDM model (e.g., 
setting Wq — —1 and Wa = 0), we determine the Hubble con- 
stant to ^ 0.025%. As recently emphasized in 144], precision 
measurements of the Hubble constant can be a key component 
of dark-energy studies; BBO would provide the most precise 
measurement of Hq that has ever been contemplated. 

In addition to the Hubble constant, BBO will directly con- 
strain the dark-energy equation of state. Figure [8] shows the 
BBO constraint on wq and Wa, for our fiducial binary sam- 
ple, with the inclusion of Planck CMB priors. We find a 
^ 0.01 constraint on wq and a ~ 0.1 constraint on Wa- We 
note that we have not assumed a flat Universe in these fits, 
nor do we incorporate any other cosmological measurements 
(beyond Planck). For comparison, we consider the stage IV 
dark-energy missions (supernovae, baryon acoustic oscilla- 
tions, and weak lensing), as listed by the dark-energy task 
force f45] , representing the state of the art in future dark- 
energy missions. The combination of all stage IV missions 
improves the task-force figure of merit by a factor 8 to 15 
with respect to stage II missions (see pp. 18-20 and pp. 77- 
78 of ll45l ). For comparison, BBO finds an equivalent figure 
of merit enhancement of ^ 100, roughly an order of magni- 
tude better than all of the stage IV missions, combined. It is 
also to be emphasized that there are still fundamental concerns 
regarding possible systematic errors in all of the stage IV mis- 
sions, and thus their combined figure of merit is undoubtedly 
optimistic. As discussed above, we expect the systematic er- 
rors associated with BBO measurements to be negligible, as 
it should be possible to build BBO such that calibration errors 
are much smaller than 10^*. 



B. Weak Gravitational Lensing and Growth of Structure 

In addition to providing precision measurements of the fun- 
damental cosmological parameters (Hq, ri,„, fi^, wq, and 
Wa), BBO will also directly measure the effects of gravita- 
tional lensing, and thus place strong constraints on the pri- 
mordial dark matter power spectrum, P{k), and the growth 
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FIG. 8: Top: Measurement accuracy of the Hubble constant, h, and 
the dark-matter density, Sim- The sohd and dashed curves map the 
1(7 and 2(7 contours, respectively. The red star denotes the true un- 
derlying model. Bottom: Measurement accuracy of the dark-energy 
equation-of-state parameters wq and Wa- 



of structure. The growth of inhomogeneities is particularly 
sensitive to gravity, and thus is a powerful way to constrain 
theories that modify gravity as an alternative to assuming a 
dark-energy component. 

One of the most powerful ways to measure the growth of 
density perturbations is through gravitational lensing shear 
maps. This is done by observing the shapes of large numbers 
10^) of background galaxies, and measuring the subtle cor- 
relations in the shapes of these galaxies due to the shear from 
gravitational lensing. The shear power spectrum at any red- 
shift is sensitive not only to the distances between observer, 
lens, and source (and thus, to the dark energy component), but 
also to the distribution of lenses. This lens distribution is a di- 
rect measure of the dark matter power spectrum as a function 
of redshift, which is in turn sensitive to the growth function of 
perturbations, and thus the gravitational force. 

BBO would provide definitive measurements of the gravi- 
tational lensing convergence power spectrum, comparable to 
state-of-the-art proposed measurements of the lensing shear 
power spectrum. BBO measures an absolute luminosity dis- 
tance to each of the ^ 10'^ binaries. The error on this mea- 
surement is almost entirely dominated by the effects of grav- 
itational lensing magnification. Once the average luminosity 



distance-redshift relation is determined (as discussed in the 
previous section), it is possible to measure the deviations from 
the background relation. Because the intrinsic uncertainty 
in the distance measured by BBO is negligible when com- 
pared with lensing (see Fig. [5]l, each individual binary thus 
becomes a direct measure of the gravitational lensing magni- 
fication along the given line of sight. The population of bina- 
ries thus provides a few times 10^ individual measurements of 
the magnification out to z 3. By evaluating the two-point 
correlation function of these magnification measurements, it 
is possible to directly measure the convergence power spec- 
trum (which is equivalent to the shear power spectrum; con- 
vergence, K, is related to magnification, fi, hy ji — 1 + 2k in 
the weak lensing limit). This approach has been discussed, for 
the case of Type la supernova distance measurements, in f46l. 
Here we follow an identical approach, using binary standard 
sirens instead of supernova standard candles. In our case each 
individual distance measurement is at least an order of magni- 
tude better, and we have an order of magnitude more sources, 
even compared to the very ambitious supernova sample con- 
sidered in |46|. We note that in what follows we focus on 
the weak lensing power spectrum, and for simplicity neglect 
strong lensing. The latter will be discussed in more detail in 
Section HVCl 

In the Introduction we provided a rough estimate that BBO 
could measure weak lensing (WL) with SNR of ~ 2 x 10^ 
for its NS-NS dataset and also - 2 x 10^ for its BH-BH 
dataset, for a total SNR of - 3 x 10^. The IDEM design 
has not yet been determined, and the WL capability of the 
mission varies quite significantly over the range of possibili- 
ties. The designs that are best-suited for WL measurements 
contain ~ 5-6 x 10^ pixels in the focal plane and would have 
a goal of measuring galaxy ellipticities to ~ 0.1%, and thus 
would require ellipticity coiTelation measurements on ^ 100 
galaxies to measure the WL effect to SNR of order 1 . (This is 
because galaxies typically have intrinsic ellipticites e ^ 0.3, 
while the correlated ellipticity due to WL is a factor ^ 10 
smaller, and SNR builds up as the square root of the number 
of galaxies observed.) Ideally, IDEM would measure shear 
for ^ 10^ galaxies, covering ^ 10* deg^ on the sky, lead- 
ing to a total SNR of psNo ~ 3000. We note that LSST is 
expected to measure weak lensing for ~ 2 x 10^ galaxies, 
out to z = 3, over ^ 2 x 10* deg^, and is thus comparable 
to the most optimistic space-based lensing missions. These 
estimates of the power of weak lensing shear measurements 
assume that systematic errors (including telescope distortion, 
shear calibration, point-spread-function correction, and red- 
shift calibration) can be beaten down to the ~ 0.1%, which is 
quite optimistic (and far better than is cuiTently possible) |47|. 

The two methods of measuring WL are rather different — 
individual magnification measurements versus correlated el- 
lipticity measurements — and a proper Fisher-matrix calcula- 
tion is required to accurately compare the science yield from 
either method. Such a calculation for BBO is now underway 
and will be published in a follow-up paper. But, crudely, we 
expect the ratio of cosmological parameter estimation errors 
to be comparable to the ratio of SNRs for the two methods, 
which is of order one, when BBO is compared to IDEM mis- 
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sions designed to maximize the WL science. 

We next calculate how accurately BBO could measurement 
the convergence power spectrum. We following closely the 
approach of |46|. The weak lensing angular power spectrum 
for magnification can be written as 



W'^(r) I 



"A 



where 



W{r) ^ 3 / dr'n{r')nr. 



Hi 



dA{r)dA{r' - r) 



c^a{r) da{r') 



(17) 



(18) 



Here r is the comoving distance, dA is the angular diame- 
ter distance, n{r) is the number density of binary systems 
(normalized so that J drn{r) — 1), and Pdm is the three- 
dimension dark matter power spectrum (calculated following 
the approach of |48 |). The error on the measurement of the 
magnification power spectrum is given by 



(2/ + l)/,kyAZ 



^binar 11 



(19) 



where /g^y is the fraction of the sky covered by the survey, 
is the binning width in multipole space, is the RMS 
uncertainty of the magnification measurement from each bi- 
nary, and A^binarics IS the surface density of the binaries. The 
first term represents the error from cosmic variance, while the 
second term represents the error from shot-noise. 

Since the BH-BH merger rate is poorly known, to be con- 
servative in the remainder of this section we shall consider 
WL measurements of the NS-NS population only. (Similar 
calculations for the BH-BH case, for a range of rates, will be 
published in later work.) In Fig.|9]we show BBO's projected 
measurement of the weak-lensing magnification map for NS- 
NS mergers. We note that the error bars are for each individual 
I mode; no binning has been performed in this figure. Defin- 
ing the signal-to-noise ratio for this measurement as 
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(20) 



we find S/N w 120, over an order of magnitude improvement 
over the equivalent measurement using 10,000 SNe over 10 
deg2 146 1 . 

We emphasize that our technique for measuring WL 
through the magnification of GW sources is entirely inde- 
pendent from the more traditional WL shear measurements, 
which observe the correlations in galaxy shapes. Since the 
systematics are a major source of concern for shear measure- 
ments, an equally powerful but fundamentally independent 
technique for measuring WL would be highly desirable. It 
will also be of interest to cross-correlate the two independent 
measures of WL (amplification and shear). This will enable 
interesting consistency tests, and directly test for the presence 
of systematics. In addition, fundamental tests may be pos- 
sible, since the relation between shear and amplification is a 
robust, and thus far untested, prediction of general relativity. 




FIG. 9: BBO's measurement of the gravitational lensing convergence 
power spectrum, based just on NS-NS mergers. The red curves show 
the error bars, for each individual / mode. Each binary is a measure of 
the gravitational-lensing magnification (convolved with the intrinsic 
error) along that given line of sight. By observing many binaries, 
BBO produces a "magnification map" of the sky. This plot shows 
the power spectrum of this magnification map, which is sensitive to 
the growth of inhomogeneities in the Universe, as well as potential 
modifications to gravity. 



C. Calibration issues and galaxy misidentifications 

We have argued that BBO has revolutionary potential as 
both a dark energy and a gravitational lensing mission. There 
are a number of potentially important systematics which must 
be addressed, foremost of which are ensuring calibration ac- 
curacy and the potential misidentification of host galaxies 
(and hence redshift of the binaries). We comment on these 
in turn. 

The analysis in this paper is based on the basic BBO design 
put forth in the BBO Concept Study 1 1 1. That design was ex- 
tremely "LlSA-like" in that the test masses are freely floating, 
with no forces applied along the arm directions. In a later pa- 
per, Harry et al. B9l pointed out a flaw in the Concept Study 
design — the laser power arriving at the photodiodes would 
saturate them — and proposed a shift to a more "LlGO-like" 
design, with forces applied to the test masses parallel to the 
arm axes, to keep the photodiode operating near a dark fringe. 
Because it is difficult to measure this applied force accurately, 
this redesign would compromise the self-calibrating quality 
that is one of LISA'S strong points. LlGO's strain calibration 
is accurate to within ~ 8%, which is well below the desired 
level for the standard siren measurements considered here. An 
author of [49 1 has indicated that calibration issues were not a 
major consideration during the re-design; until the work de- 
scribed in this paper, the motivation for a highly accurate BBO 
calibration has not been recognized |50|. We have consulted 
with interferometry experts, and they suggested several plau- 
sible solutions to the saturation problem that would preserve 
the LlSA-like calibration accuracy |51|. One possibility is to 
use optics to widen the beam, spreading the interfered light 
over an array of photodiodes. Another possibility is to keep 
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the latest LIGO-like design, but introduce an additional small 
cavity behind the test mass to calibrate the force applied ISD . 
Our confidence that a good technical solution can be found 
is enhanced by the fact that LlSA-level calibration accuracy 
is actually overkill for BBO; relative errors of 10^'', and per- 
haps larger, are acceptable for the cosmological applications 
described here. 



We also note that the current Decigo design is decidedly 
"LIGO-like" (or "TAMA-Uke"), with far greater laser power 
in the arms than for BBO, and so would have to be modified to 
achieve the calibration accuracy required for doing ultra-high 
precision cosmology. We hope that the prospect of precision 
cosmology will spur instrumentalists to publish improved de- 
signs for BBO and Decigo that have (or approach) LlSA-level 
calibration accuracy. 

Another possible source of error arises from the misiden- 
tifications of the galaxies hosting the compact-binary merg- 
ers. Misidentifications could arise for several reasons: the 
field is particularly crowded; the position error is several a 
(for 3 X 10^ sources, one expects position errors to range up 
to ^ 4.5(t); or the host galaxy is very dim, and so the bi- 
nary is incorrectly ascribed to a brighter galaxy within the 
error box. There are strategies for mitigating some of these 
effects. For instance, misidentifications would naively seem 
to be most troublesome when they lead to large outliers, but, 
for estimating cosmological parameters like Ho, wq, and Wa, 
one could clearly employ robust statistics that diminish the 
effect of such outliers. Also, since for each detection one 
will know the size and shape of the error box and the den- 
sity of galaxies within it, one should be able to re-weight the 
data based upon the confidence of the identification. This re- 
weighting would enhance the importance of BH-BH mergers 
and gamma-ray bursts, since in those cases the error volumes 
will be extremely small. Of course, misidentifications will 
lead to some degradation of BBO's WL results as well, but 
the impact here is probably less, since a good fraction of the 
WL SNR will likely come from BH-BH mergers, and BBO's 
angular resolution for these is ~ 25 times better (in solid an- 
gle) than for the NS-NS case. 

The "dim-galaxy" problem is potentially more trouble- 
some, since one is more likely to miss a galaxy on the "far 
side" rather than the near side of the error volume, which 
could lead to bias. One could attempt to quantify this effect 
by comparing results for varying exposure times (since in the 
limit of infinite exposure time, the dim-galaxy problem goes 
away). We suspect that this bias is quite small, and it is to be 
emphasized that this is a bias on the host identification (and 
hence redshift determination), not on the primary observation 
of the GW binary. Addressing this potential bias requires the 
development of a fairly detailed plan for searching for optical 
counterparts, plus a detailed, robust data analysis algorithm 
that mitigates misidentifications, which are beyond the scope 
of this paper. 



IV. FURTHER ASTRONOMY FROM BBO 

Although BBO has been primarily conceived as a detector 
of inflation-generated GWs, in this paper we argue that BBO 
would also be an unrivaled dark-energy mission. In this sec- 
tion we briefly discuss some other unique astronomical oppor- 
turuties that would be afforded by BBO. 

A. Gamma-ray bursts 

Short/hard gamma-ray bursts are widely believed to result 
from (some subset of) NS-NS or BH-NS mergers f52l. If 
this is the case, then BBO will serve as an "early warning 
system" for short/hard bursts, predicting the precise time and 
sky location of every burst, months in advance. This advance 
warning allows the bursts to be monitored with a full panoply 
of telescopes before they burst, and will permit searches for 
any "pre-burst" electromagnetic activity. BBO will also tell 
us very precisely the masses of the two bodies and the geome- 
try of the system, which one will be able to correlate with the 
electromagnetic signals. This should fully resolve the ques- 
tion of short/hard gamma-ray burst progenitors, as well as elu- 
cidate other important issues (such as the nature and amount 
of beaming). In fact, BBO may potentially overpredict the 
short/hard bursts, since not all mergers will necessarily lead 
to observable bursts. Our knowledge of the orbital geometry 
of the inspiraling binary should help predict which mergers 
will be observable (e.g., if the gamma-ray beaming is perpen- 
dicular to the orbital plane). BBO would allow us to com- 
pletely characterize the properties of orphan afterglows, re- 
sulting from GRB events beamed such that they are not visible 
to us. In addition, BBO would tell us which compact-binary 
mergers do not lead to observable bursts and/or afterglows. 

It is to be noted that gamma-ray bursts will also serve as 
"verification sources" that confirm that BBO is working as 
expected: short/hard bursts should all "go off" with merger 
times and sky-locations that are consistent with BBO's very 
precise predictions. 

B. IMBHs to high redshift 

Another interesting (but more speculative) BBO source is 
the merger of intermediate-mass black holes (IMBHs) at high 
redshift. It is predicted that the death of Population 111 stars 
could result in BHs weighing a few hundred solar masses |53 1. 
Some of these IMBHs may have grown (presumably mostly 
by gas accretion) into the very massive BHs that now exist in 
the nuclei of nearly all large galaxies. In the early build-up of 
galaxies from smaller subhalos, when the subhalos merge the 
IMBHs they contain may be expected to merge as well. BBO 
would detect these BH mergers with very high SNR, out to 
z = 20 and beyond. In Table|Il]we list BBO's SNR for IMBH 
mergers at 2 = 20, for several mass combinations. 

Note that for binaries with GW frequency at the last stable 
orbit satisfying Jlso ?i 0.5 mHz (i.e., lying rightwards of 
the minimum in the BBO noise curve), one expects the SNR 
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Ml 


le2 


3e2 


3e2 


leS 


le3 


le3 


Ma 


le2 


le2 


3e2 


le2 


3e2 


le3 


med SNR 


1.4e3 


1.9e3 


2.7e3 


1.6e3 


2.2e3 


2.2e3 



TABLE II: Median matched-filtering SNRs for inspiralling 
intermediate-mass black hole binaries (IMBHs) at redshift z = 
20. The masses are the locally measured ones (i.e., not redshifted 
masses), given in units of Mq. 



be accurate to '--^ 0.1 sec, a fractional accuracy of better than 
one part in ~ 10*. Thus independent estimates of the Hubble 
constant from time delays will be possible for ~ 300 multiply- 
imaged mergers, modulo the standard difficulty of accurately 
modeling the density profiles of the lensing galaxies. In addi- 
tion, the relative magnifications should help constrain the dark 
matter density profiles, as well as directly break the mass- 
sheet degeneracy. 



of the inspiral waveform to scale like . However, requir- 
ing !lso > 0.5 mHz requires M(l + z) < 2.2 x lO^M©, 
which is satisfied only by the first column in the table. That is, 
the tabulated SNRs are the results of two competing effects: 
as the masses increase, so does the signal amplitude, but si- 
multaneously less and less of the inspiral waveform lies in the 
sensitive portion of the BBO band. 

M. Volonteri has supplied us with the results of her merger- 
tree model, based on two assumptions: 1) today's massive and 
supermassive BHs began as initial seed BHs of a few hundred 
Mq, and 2) these seeds formed at z > 20, and the seed BHs 
grew efficiently from accretion. Using Volonteri's model, we 
estimate that BBO's detection rate for IMBH mergers in the 
mass range Mi + M2 < 1,OOOM0 would be 30/yr, of 
which ^ 25/yr would be at z > 10. By comparison, Sesana 
et al 1 54 1 calculate that a network of three third-generation 
ground-based interferometers would detect IMBH mergers at 
a rate of ^ 2/yr (based on a very similar merger-tree model 
of Volonteii), and almost all of these detected IMBHs would 
be at 2 < 8. Thus BBO offers a unique opportunity to directly 
observer the very first seed black holes in the universe. 



C. Strong Lensing 

In addition to a direct measurement of the convergence 
power spectrum, BBO will measure the full magnification 
probability distribution function as a function of redshift, in- 
cluding resolving the high-magnification tail. These high- 
magnification effects can be very sensitive to the dark-matter 
(and baryon) profiles in galaxies (see, e.g., ||55l for an optical 
version of this). For example, galaxies with more concen- 
trated, cuspy profiles may be more likely to engender strong 
lensing than those with more flattened cores. 

The fraction of multiply-imaged quasars is ^ 10^"^. One 
expects a comparable fraction of multiply-lensed NS-NS bi- 
naries: i.e., ^ 300 multiply-imaged binaries. In the BBO 
dataset, these will appear as pairs of binaries having nearly 
identical sky locations and essentially identical masses, spins, 
and orientations (i.e., identical to within the error bars), but 
with different apparent distances (according to the magnifi- 
cation of each image), and with arrival times differences of 
order months to years. The observed rate of multiple-imaging 
is an important probe of dark-matter density profiles, as well 
as the overall dark matter distribution (e.g., erg), and BBO will 
provide an extremely clean measurement. BBO's determina- 
tion of the time-delays between the multiple images should 



V. SUMMARY, CONCLUSIONS AND FUTURE WORK 

Studies of the dark energy generally take one of two ap- 
proaches: The first approach is to measure the luminos- 
ity distance-redshift relation (or angular diameter distance- 
redshift) to high accuracy. Type la supernovae and baryon- 
acoustic-oscillation measurements fall into this category. The 
second approach to exploring the dark energy is to measure 
the weak-lensing shear power spectrum, and infer values for 
the growth of structure and cosmological distance ratios. In 
this paper we have shown that BBO will provide unprece- 
dented measurements of both the luminosity distance-redshift 
curve as well as the weak lensing power spectrum. 

As discussed in Secjlj the success of BBO as a cosmolog- 
ical probe is dependent upon obtaining redshifts for a large 
number of the host galaxies to the binary sources. It is to 
be emphasized that any subsample of galaxies is sufficient 
(e.g., LRGs), since there are no standard siren systematics 
expected to be associated with the nature of the host galaxy. 
Although obtaining redshifts for 3x10^ host galaxies is cer- 
tainly a demanding requirement, as mentioned in Sec. |lj we 
expect that surveys such as LSST or BigBoss may provide 
the required dataset in due course. Future work will be re- 
quired to refine the estimates of BBO's performance given 
here, and to consider ways that the mission design might be 
modified to improve its price/performance ratio from a cos- 
mological perspective. We hope that this work will encour- 
age GW instrumentalists to develop workable designs with 
calibration accuracy comparable to LISA'S. We plan on im- 
proving the analyses presented here in several ways. First, in 
this paper we made simplifying approximations in calculat- 
ing BBO's parameter estimation accuracy. Although we ex- 
pect that these simplifications will have little impact on the 
results, we will explicitly check this by 1) including the effects 
of non-zero eccentricity and Lense-Thirring precession in our 
waveform model, 2) going beyond the Fisher matrix approx- 
imation in calculating the expected sizes of the eiTors |56|, 
and 3) calculating how cross-coiTelations between different 
binaries affects parameter estimation for each. Regarding the 
cosmological constraints derivable from BBO's dataset, we 
intend to 1) explore BBO's performance as measured by other 
dark-energy figures of merit, 2) improve the sophistication of 
our method for incorporating priors (and, in particular Planck 
CMB priors) into the BBO analysis, 3) incorporate robust pa- 
rameter estimation methods that mitigate the effects of host 
galaxy misidentifications, and 4) generalize our approach to a 
model-independent, multi-parameter description of the dark 
energy ll57l l58l |59l , instead of the arbitrary two-parameter 
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form (wq and Wa) considered here. In addition, we will inves- 
tigate whether there are special synergies in cross-correlating 
BBO's gravitational lensing magnification measurements with 
the weak-lensing shear measurements made by a JDEM-like 
mission. We also plan to further explore BBO's potential 
impact upon other areas of astronomy (some of which were 
sketched out in Sec. |lV| i. Lastly, we will investigate how var- 
ious figures of merit vary with BBO's sensitivity and other 
mission design parameters. There are simple variants to the 
current design which may retain many of the high-precision 
cosmology applications while considerably reducing the price 
or risk of the mission. For instance, having four constellations 
of mini-LlSAs is overkill for high-precision cosmology; three 
constellations are certainly sufficient, and two may be suffi- 
cient, if combined with measurements from planned, third- 
generation ground-based GW detectors (such as the Einstein 
Telescope ll60l ). In the latter case, the mission would require 
six satellites instead of twelve. And while a noise curve at 
the level shown in Fig.|2]is probably necessary for detecting 
a stochastic GW background, it is somewhat overqualified for 
the task of doing ultra-high precision cosmology: as currently 
conceived BBO could detect all NS-NS mergers out to z = 5, 
while for doing high-precision cosmology it is probably suf- 
ficient to detect and localize a reasonable fraction of mergers 
out to 2; '--^ 2. Since BBO would surely be a multi-billion dol- 
lar mission, it is important to look for ways of significantly 
reducing mission cost while retaining a large fraction of the 
science. 

In this paper we have shown that BBO is a particularly 
powerful mission, as it will provide revolutionary measure- 



ments of both the luminosity distance-redshift relation and 
the growth of structure (through gravitational lensing mea- 
surements). Unfortunately, a deci-Hz GW mission like BBO 
is, optimistically, at least fifteen years from being built. Nev- 
ertheless, the power of a BBO-like mission for precision 
cosmology — stemming from very precise, unbiased distance 
measurements of '-^ 3 x 10^ NS-NS binaries to z '--^ 5 — is 
so revolutionary that BBO could represent the future of high- 
precision cosmology. 
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